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New criteria for the equation of state development: 
Simple model fluids 

Yurko Duda and Pedro Orea 1 
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Eje Central 152, 07730 Mexico D.F., Mexico; e-mail: yduda@lycos.com 

Recently we have proposed (J. Chem. Phys. 128 (2008) 134508) a new rescaling of fluid density 
p by its critical value pj to apply the corresponding states law for the attractive Yukawa fluids 
study. Analysis of precise simulation results allows us to generalize this concept to the case of simple 
fluids with different interparticle interactions, like Mie (n,m) and Sutherland pair potentials. It is 
shown, that there is a linear relationship between the critical pressure and critical temperature, as 
well as the critical density and inverse critical temperature for these frequently used pair potentials. 
As a consequence, the critical compressibility factor of these model fluids is close to its universal 
value measured experimentally for different real substances. 
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Equation of state (EOS) is of crucial importance 
to predict or correlate fluid thermodynamic properties 
[l|, y, |3j . Since van der Waals proposed the first version 
of his EOS [J], many modifications have been proposed, 
and today, there are numerous more accurate EOSs re- 
ported in the literature for describing the behavior of 
simple fluids and fluid mixtures. An improvement in the 
accuracy of analytical EOSs for the modeling of exper- 
imental data of real fluids is usually achieved, not by 
improving theoretical base, but by introducing empirical 
temperature functionalities in their parameters. Obvi- 
ously, such empirical or semi-empirical approaches, be- 
side to require a lot of experimental [l|, |2| or simula- 
tion data [a, |6(, forbidden molecular interpretation of 
results. On the other hand, theoretically based EOSs 
are very general and applicable for the entire family of 
systems obeying t he g iven interparticle potential model 
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The most apparent progress toward EOS was made 
by applying pri ncip les of statistical mechanics. In- 
tegral equation |i p| . 11 1, thermodynamic perturbation 
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161 ], and mean field theories 
the most common theoretical approaches to develop an- 
alytical EOS for simple fluids. The advantage of the- 
oretically based EOS over their empirical counterparts 
is that they can be further improved by testing against 
computer simulations and later be used as a reference 
for perturbation theory calculations of bulk and intcrfa- 
cial properties of more complex fluids such as polymers 
and associating systems (i.e., hydrogen-bonding fluids) 

[S El El El, m, Hi. 

Although last decades significant progress has been 
achieved in the development of theoretically based ex- 
pressions for EOS, there is still a lack of equations that 
described PVT properties of simple fluids in wide range 
of thermodynamic parameters. Such precise analytical 
expression may be obtained if one has a set of univer- 
sal criteria (as for example, corresponding states law 
[24 I25L |26J ]) or/and apply self-consistent integral equa- 



tions [27] . 

Recently, we have applied the corresponding states 
law for the hard core attractive Yukawa (HAY) fluid 
[26j, and revealed that after a certain rescaling of fluid 
number density p the reduced surface tension and pres- 
sure map onto the master curves. To make such anal- 
ysis we performed canonical Monte Carlo (MC) simu- 
lations with relatively high precision [28|. The goal of 
the present work is to show that the master curve for 
the supercritical reduced pressure is not restricted to the 
HAY fluid and is universal for other simple nonelectro- 
static fluid pair potentials. Namely, in this communica- 
tion we consider the following widely used pair potentials 

UlilHEl,!!!,!!!: 



1) Mie (n,m) potential, 



u(r) = e 
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2) Sutherland potential, 
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if r > 1, 
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and, 



3) HAY potential, 



u(r) 



00, if r < 1, 

— eexp[— «(r — l)]/r, ifr>l, 



(3) 
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where e is the potential well depth, k defines the range 
of the potential, and r is an interparticle distance re- 
duced by particle diameter, which is chosen to be a unit 
of length. 

Note, that last two potentials, eqs. flU) and ([3]), have 
hard sphere repulsive part, while the family of Mie po- 
tentials (which includes widespread used Lennard-Jones 
(LJ) pair potential, when n = 12, and m = 6) reproduces 
a soft repulsive part at short distances and a smooth dis- 
persion part at intermediate and long distances. In other 
words, the distinct nature of the potentials described 
above is obvious. 

In this work we mainly use the simulation data of crit- 
ical parameter values reported in the literature 28|, 129 . 
30 . l3lL |32j • Some new critical data shown in Table 1 were 
estimated by fitting the MC simulation data for coexist- 
ing liquid pl and vapor py densities 34] to the rectilinear 
diameter law 33J, [35j 



PL + Pv 



Pc + A\T/T c -l\, 



and to the following equation 



PL-pv = B\T/T c -l\ 



(4) 



(5) 



where A and B are the coefficients of the leading terms 
in the general Wegner expansions, and f3 — 0.325 is the 
universal value of critical exponent in the case of short- 
ranged potentials 1331 pal , that belong to the Ising uni- 
versality class [3a, |37|. Up to four pairs of densities, 
Ph{Ti) and p v (Ti), at temperatures T, close to its criti- 
cal value, have been used in Eqs. dU) and (0 to estimate 
the critical parameters p c and T c ; averaging of the six 
values (combination of the four pairs of density) has pro- 
vided us the mean values and uncertainties of the critical 
values presented in Table 1. It is well accepted [36] that 
such kind of estimation gives the same values of critical 
data as obtained by other more sophisticated simulation 
techniques [33ll37|. 

The critical pressures were estimated on the base of the 
Clasius-Clapeyron equation. Besides, we have performed 
canonical MC simulations to calculate supercritical pres- 
sures for the fluids with pair potentials defined above. All 
details of the simulation can be found elsewhere 261 12a] ; 
as usual, additional simulation runs of different size sys- 
tems and/or with different length of r cu t have been per- 
formed to verify that our results do not depend on these 
parameters. 

First of all, let us consider the relation between the 
critical parameters of the liquid-vapor phase transition 
of the fluids with different pair potentials. As seen in 
Fig.l, the critical pressure P c is a linear function of the 
critical temperature T c for the different potentials consid- 
ered. Note that family of Mie potentials is represented 
by more than 20 points which correspond to the different 



combinations of parameters n and m in eq. iJTJ . Besides, 
in Fig.l there are some points of the slightly modified Mie 
potentials, which were studied in Ref. [32j; namely, the 
Buckingham exponential-6 (Bexp6), and the LJ(a — 6) 
potentials. It is important to remark, that the linear 
function P c — P C (T C ) is nearly the same as was reported 
in Ref. [26j for HAY fluid potential with different values 
of K. 

As shown in Fig. 2, the critical densities p c are linearly 
dependent on the inverse critical temperature. Some of 
the data on the p c = p c (l/T c ) curve (for example, the 
data of Ref. [32j , which were not used for the linear fit- 
ting) are located slightly out of the adjusting solid curve, 
although presenting linear dependence too. It is impor- 
tant to note, that the two linear functions shown in Figs. 
1 and 2 allow to obtain the value of the critical com- 
pressibility factor, Z c = P c /(p c T c ), around 0.30 ± 0.01 
for all the model fluid considered here. This value of Z c 
is close to the universal value reported for different real 
substances [24|, and usually is being used for the critical 
pressure estimation. Prediction of Z c obtained from var- 
ious theoretical approaches 1CJ, [l3j, [l5j, |lg, [l7|, [l8|, [l9( are 
far from this value because, as well known, none of them 
are sufficiently precise in the critical region. 

Finally, in Fig. 3 the reduced pressure, Pr — P/P c , 
as a function of the rescaled density, p r = pj ' pj , is 
presented for the four different potentials. Such rescaling 
of density has been proposed and tentatively justified in 
our previos work [261 ] , where the extended corresponding 
states law has been applied for the description of the HAY 
fluid properties. As can be seen in the figure, for the 
three reduced temperatures considered, all four model 
fluids match on the same master curves at each value 
of Tr = T/T c . The best agreement of the simulation 
data is reached at the lowest temperatures and densities; 
only slight discrepancy (almost within the error bars) is 
observed for the highest densities (or pressures). 

In summary, novel criteria for estimating the critical 
parameters and reduced equation of state of simple fluid 
model potentials is proposed. Namely, the rule of linear 
dependence of critical pressure on critical temperature, 
as well as critical density on inverse critical temperature 
is tested for a number of potentials. Such rule gives criti- 
cal compressibility factor around 0.30 ±0.01 which agrees 
well with experimental data of some real gases. Besides, 
new re-scaling of simple fluid density, p r = pj pj allows 
one to describe the reduced pressure by a single reduced 
isotherm. We believe, that the proposed universality re- 
lations will help to improve old theoretical approxima- 
tion in the field of simple fluid modeling, and develop 
new ones. Moreover, it would be interesting to analyze 
the available experimental PVT data with the idea to 
verify the applicability of the corresponding states rules 
presented here. More extensive computer simulations are 
necessary to investigate their applications for the inter- 



facial properties of simple model fluids. This issue will 
be considered in future reports from our laboratory 34 1 . 
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Figure Captions 

Fig. 1 Relation between critical pressure P c and 
temperature T c for different model potentials: HAY with 
1.5 < k < 7 0, LJ(a-6) with 10 < a < 22 ^, Mie£n-6) 
with 7 < n < 32 [H, Mie(2n - n) with 6 < n < 18 [p], and 
Bexp6 with 12 < a < 22 gl|. 

Fig. 2 Relation between critical number density p c 
and inverse critical temperature T c _1 of the model potentials 
considered in fig.l, and Sutherland potential with 7 = 4, and 
6 [1} 
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Fig. 3 Reduced pressure, Pr 
of reduced fluid density p r = pj ' pj . 
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TABLE I: Critical parameters of the four model potentials 
obtained from the MC simulations, and used for the Pr = 
Pn{pr) function calculations (see Fig. 3). The last column 
presents the cutoff distance r cut applied for each pair poten- 
tial. 

Potentials T c p c P c Z c r cut 

HAY k, = 1.8 (Ref. [28]) I.I8O5 

Mie (12-6) (LJ fluid) [34] 1.290 9 

Mie (32-6) [34] 0.865 4 

Sutherland (7 = 6.0) [34] 0.595 5 
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